DISTINCT OXYGEN HOLE DOPING IN DIFFERENT LAYERS OF 

Sr2Cu04_5 - La2Cu04 SUPERLATTICES 
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X-ray absorption in Sr2Cu04_i — La2Cu04 (SCO-LCO) superlattices shows a variable occupation 
with doping of a hole state different from holes doped for x < XopUmai in bulk La2-xSrxCu04 and 
suggests that this hole state is on apical oxygen atoms and polarized in the a — b plane. Considering 
the surface reflectivity gives a good qualitative description of the line shapes of resonant soft X-ray 
scattering. The interference between superlattice and surface reflections was used to distinguish 
between scatterers in the SCO and the LCO layers, with the two hole states maximized in different 
layers of the superlattice. 
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I. INTRODUCTION 

Hole doping of La2-xSrxCu04 is described in the 
Zhang-Rice singlet (ZRS)-upper Hubbard band (UHB) 
model for < a: < 0.2 as ZRS states on in-plane oxy- 
gen atoms. These holes are visible as a feature in X-ray 
absorption spectroscopy (XAS) at the O K edge, called 
the mobile carrier peak (MCP). The MCP gains inten- 
sity with doping at the expense of LHB (lower Hubbard 
band) states and concurrent "spectral weight transfer" 
from UHB states. This model, expanding the Cu d^'^-y'^ 
one-band Mott-Hubbard model to the Cu dx2-y2, O Px, 
O Py three bands of the Cu02 planes, describes well the 
variation, observed with XAS at the O K edge, of un- 
occupied density of states in cuprates for relatively low 
doping (x < 0.2)L2. 

However, the variation of the maximum critical tem- 
perature Tc^max between different superconducting com- 
pounds cannot be explained within the UHB-ZRS model 
and its in-plane orbitals only and possible extensions of 
the model to out-of-plane orbitals have been intensively 
investigated. Since there can be only one Fermi surface 
for an isolated Cu02 plane, other indications of the rel- 
evance of out-of-plane orbitals came from angle-resolved 
photoemission (ARPES) measurements of square- and 
diamond-shaped Fermi surfaces for La2Cu04 (Ref. 3) 
and Ca2Cu02Cl2 (Ref. 4). The most relevant out-of- 
plane orbital that hybridizes with states in the Cu02 
plane is the apical oxygen pz orbital mixed with the Cu 
3d322_^2 (Refs. 5,6) or 4s (Ref. 7) orbital in the Cu02 
planes. The occupation of the apical oxygen orbitals 
modifies the bond valence sums and was used to explain 
general trends of Tc,max'^'^ The effect of apical oxygen pz 
energy level on parameters of an expanded t — J model 
has been considered for different materialsJ^ In contrast, 
the axial hybrid between the O Pz and Cu 4s orbitals, 
with Cu (i322_r2-states occupied, has been addressed in 
Ref. 7. These calculations predicted that an empty api- 
cal oxygen pz orbital modifies the in-plane hole hopping 
parameter t' between sites along the orthorhombic axes, 
consistent with ARPES measurements, with a suppressed 
t' from the presence of unoccupied apical oxygen orbitals 



correlating with a smaller T^max'^'^'^ 

In addition, the depletion of the UHB states in bulk 
La2-xSrxCu04 at a; ~ 0.2^ shows that the ZRS-UHB 
model also needs to be modified to describe the hole dop- 
ing for higher x. However, because of limits of bulk crys- 
tal growth, doping dependence studies have been limited 
to doping near x = 2, realized in bulk Sr2Cu04_5 with a 
very large T^r^— where the oxygen atoms are removed 
from the structure (J > 0), and for relatively low x, near 
the superconducting dome, where there is evidence of 
a qualitative change with x in the doping process near 
^optimal- Specifically, the effective Cu ion magnetic mo- 
ment in La2-xSrxCu04 and the magnetic exchange be- 
tween Cu spins are strongly reduced^ at x ~ 0.2 and 
the T-independent Pauli paramagnetism is replaced at 
x ~ 0.22 by T-dependent Curie paramagnetism with in- 
creased doping. Calculations^'^ suggested that the dop- 
ing mechanism is different for x > XopUmaU as the num- 
ber of ai-symmetry states, related to the Cu Sd^z^-r^ 
orbital^° and exceeding the optimal doping XopUmai, cor- 
related with Tc,max- A small but growing contribu- 
tion from apical oxygen pz orbital holes mixed with Cu 
d^z^-r^ orbitals at MCP was inferred from angle- resolved 
XAS measurements^^ ^ and from the variation of the Cu- 
apical oxygen distance with doping2^, which also corre- 
lates with a variation in T^^^n^ Recently, a saturation 
of the XAS MCP intensity with increasing doping near 
^optimal was observed^^ and dynamical mean-field theory 
calculations^^ concluded that either additional orbitals 
become relevant in this doping range or that new model 
parameters would be needed to account for multiple- 
hole interactions. The additional orbital or band to con- 
sider in the La2-xSrxCu04 hole doping process between 
X ^ 0.2 and a; ~ 2 is not known. 

In this study, doping in the range 1 ^ x < 1.5, in- 
accessible with bulk crystal growth techniques, was ob- 
tained with Sr2Cu04_5 — La2Cu04 (SCO-LCO) super- 
lattice (SL) growth. We observed that occupations of 
two distinct oxygen hole states are gradually modified 
with doping. Using the interference of the SL refiections 
with the surface reflection, we determined that the two 
hole states are maximized in different layers. This sug- 
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FIG. 1: (Color online) (a) Measurements before and after 
ozone annealing (OA) show that all SL are superconducting. 
The inset shows a sketch of one SL superperiod. The atomic 
planes within the SCO and LCO layers are shown in Fig. 6(a). 
(b) Hard X-ray reflectivity characterization {E — 8.05 keV). 
The sharp peak at 20 ~ 14° is the substrate (002) reflection 
peak, with clsao = 12.64 AM Using this peak as a refer- 
ence marker, HXD measurements showed that the average 
ML thickness was 6.64 A, 6.63 A and 6.59 A for SL-A, SL- 
B and SL-C, respectively. The difference is consistent with 
the thicknesses of the SCO and LCO layers in each SL and 
the larger c-axis parameter of SCO films csco = 13.55 
compared to that of LCO films clco = 13.3 A. Superlat- 
tices with smaller LCO layer thickness in a superperiod (not 
shown) had less good growth. The inset shows the epitaxial 
growth, where aLSAO ~ 3.76 Ai^ (c) Characterization with 
soft X-rays near the La edge (E = 832 eV). The dotted lines 
show the evolution of the SL peaks with superperiod thick- 
ness.— Inset shows L scans at MCP and SHP energies (Sec. II 
D). The difference at L=3 between MCP and SHP amplitudes 
is also visible in Fig. 4. 



gests that the states emptied preferentially in the SCO 
layers are the additional states in the extension of the 
ZRS-UHB model to this doping range. 



II. EXPERIMENTS 
A. Superlattice structure 

The superlattice samples were grown by molecu- 
lar beam epitaxy on LaSrA104 (LSAO) substrates at 
Brookhaven National Laboratory. Three SL, called SL- 
A, SL-B and SL-C in the following, were chosen after 
atomic force microscopy (AFM) and hard X-ray diffrac- 
tion (HXD) measurements. HXD measurements [Fig. 
1(b)] were made using a Philips X' Pert diffractome- 
ter. From a Nelson-Riley fitting, the number of lay- 
ers in one superperiod for SL-A and SL-C is iV = 
7.81 ±0.1 ML and N = 9.27±0.3 ML respectively, where 
1 ML ( "molecular layer" ) is the average d-spacing of one 
SL "formula-unit" layer (half the unit cell), given by the 
(002) reflection. The SL structure was further charac- 
terized with resonant soft X-ray scattering (RSXS) mea- 
surements at beamline XIB at the National Synchrotron 
Light Source. Reflectivity measurements near the La 
edge [Fig. 1(c)] for Q = (0,0, Q^), where = 2ttL/csl 
is the scattering momentum in units of SL superperiod 
csL, show superperiods with an integer number of lay- 
ers, iV = 8 ML and TV = 9 ML for SL-A and SL-C re- 
spectively, consistent with the HXD measurements. The 
shghtly worse SL-B has N = 7.8 ML. The thickness of 
the SCO layers within a SL superperiod is Nsco ~ 2 ML 
for all samples, but somewhat larger for SL-A compared 
to SL-B or SL-C, because of larger doping (Sec. H B). 
The number of repeats was 8, 7 and 8 for SL-A, SL-B 
and SL-C, respectively. A sketch of one SL superperiod 
is shown in Fig. 1(a) (inset). 

The superconducting critical temperatures from AC 
susceptibility measurements were 25 K for SL-A, 38/39.5 
K for for SL-B and 14/36.5 K for SL-C, where the flrst 
and second values for SL-B and SL-C are for measure- 
ments before and after ozone annealing at 350 °C for 20 
min [Fig. 1(a)]. AFM images were taken on a Dimen- 
sion 3100 instrument. They showed a surface covered by 
islands approximately 150 nm wide and 2 nm high for 
SL-B (data not shown) ; SL-A or SL-C did not have these 
features. The AFM surface RMS roughnesses CTg for SL- 
A, SL-B and SL-C were 0.61 nm, 0.76 nm, and 0.43 nm, 
respectively. 

The interface roughness can be characterized with 
HXD reflectivity. The L^l reflection width [Fig. 1(b)], 
dominated by the total SL thickness and not sufflcicntly 
sensitive to small-scale roughness, was approximately the 
same for all SL. However, the L = 2 linewidth was the 
same for SL-A and SL-C and 9 % higher for SL-B. Also, 
the I{L — 2)/I{L — 1) ratio was the same for SL-A 
and SL-C and 23 % smaller for SL-B. These observa- 
tions suggest that SL-A or SL-C have a smaller interface 
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FIG. 2: (Color online) (a) Superlattice fluorescence yields dif- 
fer at the MCP and SHP energies only. The relative intensities 
of the MCP peaks in fluorescence yield make it clear that SL- 
A has more Sr than SL-B and SL-C. Since the LSAO substrate 
SHP is smaller than that of either SL, sampling more sub- 
strate FY in SL-A measurements (SL-A is thinner than SL- 
C) cannot explain the higher SHP intensity. The inset shows 
for comparison FY for bulk La2-xBaxCu04 (with x = 0.125), 
La2Cu04+5 (with S = 0.12) and NaxCa2-xCu02Cl2 (with 
X — 0.08). (b) Angle-resolved fluorescence yield on SL-A. 
The measurements have been aligned at 525 eV and 560 eV. 
Inset shows the measurement geometry. The detector angle 
29 was kept flxed at 110°. 



roughness than SL-B, consistent with the AFM and su- 
perperiod measurements. In addition, although the sur- 
face roughness Us of SL-A is larger than that of SL-C, 
their interface roughnesses ai are similar. This will be 
used to explain the difference in scattering at the O edge 
between SL-A and SL-C (Sec. Ill B). For an estimate of 
the Sr doping range we use a SL interface RMS roughness 
6 A. The maximum Sr doping is estimated with ei- 
ther a "flat-top" or a Poisson distribution as Xmax ~ 1-5 
and Xmax 1 for a SL with 2 ML and 3 ML thick SCO 
layers, respectively, outside the current possibilities of 



bulk crystal growth. Therefore, the middle nominally 
SCO layer in a 3 ML thick SCO layer is approximately 
Lao.5Sri.5Cu04 and the two nominally SCO layers in a 2 
ML thick SCO layer are approximately LaSrCu04. The 
shorthand notation "SCO" and "LCO" from the deposi- 
tion sequence will continue to be used for simplicity for 
the SL layers. Our conclusions do not depend on the ex- 
act values of SCO and LCO layer thickness or interface 
roughness. 



B. Oxygen edge fluorescence yield 

The first indication of hole doping at two distinct 
energies comes from XAS measurements with fluores- 
cence yield (FY) detection and 7r-polarized incident light, 
made at beamline XIB at the National Synchrotron Light 
Source. The sample and detector angles were 9 = 80° 
and 29 = 110° respectively, defined as shown in Fig. 
2(b) (inset). SL measurements are shown in Fig. 2(a). 
FY of bulk La2-xSrxCu04 (LSCO) has three main low- 
energy features^: the MCP of the ZRS state, the UHB 
feature and a peak called here the "second hole peak" 
(SHP). The highest peak at 536.1 eV, sometimes associ- 
ated with orbitals mixed with La states^^, is followed by 
"continuum oscillations" . There is no clearly discernible 
UHB peak in SL FY. The UHB intensity in bulk LSCO 
is negligible^ for x > 0.15; therefore, the nominally un- 
doped LCO layers in the SL are instead thoroughly doped 
with MCP holes (the occupation factor (tf^LCo) > 0.15), 
consistent with estimates based on roughness and hole 
diffusion length^^. 

The SL FY spectra in Fig. 2(a) have been aligned 
below and above the edge. The alignment normalizes 
out the variation in the number of oxygen atoms and 
therefore, the remaining differences between FY spectra 
are due to electronic contrast from the difference in va- 
lences of oxygen atoms. For SL-A and SL-C this occurs 
exclusively at MCP and SHP energies. Therefore, the 
SHP peak is related to doping of oxygen states, not to 
a structural defect, e. g. a vacancy or interstitial oxy- 
gen. Further support for this interpretation comes from 
FY measurements of bulk LSAO, where the SHP peak 
is observed (Fig. 2); since LSAO has no vacancies or in- 
terstitial atoms, SHP cannot be related to these defects. 
This is also consistent with the observation that SL-A 
has a larger SHP intensity than SL-C, while compounds 
without apical oxygen atoms, e. g. NCCOC, have a 
smaller intensity at SHP [Fig. 2(a), inset]. These FY 
results are a clear indication that the ZRS-UHB model 
extension in our doping range beyond the MCP and UHB 
levels is related to SHP and not to the vacancies of bulk 
SCO. 

The SHP energy is relatively insensitive to the type 
of oxygen neighbors. For instance, it is not influenced 
by the Cu d states since Al^"*" in LSAO has no d states 
or by replacing La with Sr. The FY increase from SL- 
C to SL-A is at the same energy (SHP) as the energy 
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FIG. 3: (Color online) (a) Reflectivity at the La edge for SL- 
A. The plots have been shifted vertically for clarity. Inset: 
FY for SL-A and SL-C. The only clear difference in FY is an 
increase in the overall height, consistent with a larger num- 
ber of La atoms in SL-C. The two lower curves are the LSAO 
substrate FY and reflectivity a.t 6 = 60°, showing the strong 
resonant increase in the substrate reflectivity. The difference 
from measurements on SL is due to the absence of SL scat- 
tering contrast SfsL (Sec. Ill A), (b) Line shapes at the Cu 
edge. For intermediate conditions the scattering line shape 
is flat: this is observed for SL-C at L—i (not shown). The 
two-dimensional proflle has a saddle point for L=l: it is a 
peak when scanning L (inset) and a dip when scanning E. 
This is in contrast to multiple-scattering effects in HXD that 
can be observed for certain oxide SL, where a dip is present 
at integer L instead of a peak (data not shown). The lower 
curves compare the FY for SL-A (full symbols) and SL-C 
(open circles). 



of a peak of the undoped LCO compound, which sug- 
gests that the SHP state, partly empty in undoped LCO, 
is gradually emptied further, rather then removed, with 
additional hole doping in SL-A compared to SL-C; that 
is, doping also makes holes at SHP. Therefore, the to- 
tal numbers of MCP and SHP holes in one superpe- 
riod are related by the charge conservation equation as 
J2i{tf^^) + J2i{tf"^) = 2iVsco, where I sums over 
one superperiod, A^sco is the number of SCO layers and 
2^^300 is the total number of doped holes in one super- 
period (each Sr2Cu04_5 layer dopes 2 holes when 6 ~ 0). 

Knowledge of the SCO oxygen vacancy site would sug- 
gest the site of SHP; however, the question of the vacancy 
site is not settled. The absence of this feature in the 
cuprate NCCOC, which does not have the apical oxygen 
(Fig. 2), suggests that SHP is an apical oxygen state. It 
is possible that SHP is absent in NCCOC for other rea- 
sons than the absence of the apical oxygen. For instance, 
in an alternative view of angle-resolved FY measurements 
on bulk LNO^i^, SHP are considered in-plane oxygen 
states mixed with Ni, and polarized in-plane. However, 
the SHP energy is not influenced by the Cu d states or 
replacing La with Sr, which suggests various hybridiza- 
tions are not essential. The incipient out-of-plane doping 
for low X (Refs. 2, 22) supports the view that doping 
is not confined to the Cu02 planes for high x. Indeed, 
more detailed early calculations of the La2-xSrxCu04 
system^ suggested that the shift with doping in the rel- 
ative alignment of the in-plane [0(1)] oxygen and apical 
[0(2)] oxygen ionization potentials will eventually lead 
to the preferential emptying of the apical oxygen orbitals 
at larger x. If SHP is an apical state, the Cu ions would 
be surrounded by holes on all six neighboring oxygen 
atoms. This transition to a gradual emptying of apical 
oxygen orbitals, balancing the valences of all six oxygen 
atoms near a Cu site (this hole distribution is seen in 
undoped La2Ni04i^) is a sensible intermediate step to- 
ward the process of vacancy creation in SCO. This kind of 
change with doping in the site of the doped holes has been 
seen in other cuprates, for instance in YBa2Cu306+5, 
where the holes stay in chains in Cu and O orbitals up 
to 5 ~ 0.25, only then doping the oxygen states in the 
Cu02 planes^i^i 

To address the question of the orientation of the state 
behind SHP, angle-resolved FY measurements have been 
made. A quantitative analysis of the angular dependence 
of FY requires measurements on bulk crystals cut at a se- 
ries of angles with respect to the crystallographic planes, 
to account for footprint and self- absorption effects'^". Su- 
perlattice samples cannot be grown at arbitrary angles. 
However, the SL footprint effects are almost identical for 
MCP and SHP energies because of the very similar scat- 
tering geometry. Self-absorption effects would have to be 
very strong for the observed suppression of intensity at 
SHP. That SL self-absorption effects are relatively small 
is supported by the very small difference in the momen- 
tum linewidth between MCP and SHP [Fig. 1 (c) , inset] , 
which shows that the entire SL is probed at both ener- 
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gies. In addition, the O edge step height is approximately 
the same before the ahgnment in Fig. 2(b), which also 
suggests that self-absorption effects are relatively small 
at the O edge. 

The MCP and SHP intensities follow a similar angu- 
lar dependence [Fig. 2(b)]. We assume that the SHP 
state can be associated with a single orbital (if localized) 
or a collection of orbitals of the same type (if delocal- 
ized, as for the ZRS state, corresponding to MCP), in 
which case the SHP state would correspond to a band of 
a specific symmetry. The similar variation with angle of 
MCP and SHP intensity in angle-resolved FY measure- 
ments [Fig. 2(b)] suggests that the state corresponding 
to the SHP energy is oriented as the in-plane ZRS. Also, 
the MCP and SHP scattering amplitudes remain simi- 
lar over a wide angular range (Fig. 4), which would be 
difficult to explain if they were polarized in different di- 
rections. Therefore, our data provide evidence that SHP 
is a state at the apical oxygen site polarized in the a — b 
plane. There are three relevant apical oxygen orbitals: 
Pz pointing to the Cu in the neighboring Cu02 plane, 
Px,y pointing to Sr or La in the same LaO plane and Pz 
pointing to the Sr or La in the neighboring LaO plane. 
The only possibility consistent with this interpretation is 
the orbital pointing to the La atoms in the same LaO 
plane. 



C. Fluorescence yield and scattering at the La and 
Cu edges 

RSXS can measure bulk^"-^ and SL^'^'^ charge or- 
der. SL scattering measurements were made with tt- 
polarized light at beamline XIB in an UHV diffractome- 
ter. Before the more complex SL reflectivity at the O 
edge, we present the measurements at the La and Cu 
edges, which will be used to illustrate the model of Sec. 
HI A. 

FY measurements at the La M5 edge for SL-A and 
SL-C are shown in Fig. 3(a). The scattering con- 
trast between the SL layers from L = 1 to L = 7, 
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difference between the number of La atoms in the LCO 
and SCO layers, {tfico) ^'^'^ {^fsco)- Consistent with 
this observation, there is a large non-resonant tail in scat- 
tering at the La edge for SL peaks, which was used to 
characterize the SL structure [Fig. 1(c)]. The line shape 
in Fig. 3(a) is approximately the same for different L 
because of the relatively small contribution from the sur- 
face reflectivity compared to the SL reflection (Sec. HI 
A). 

FY measurements at the Cu edge [Fig. 3(b)] are very 
similar for SL-A and SL-C; the scattering contrast at 
this edge is given by a small difference in the Cu valence, 

= Usco - hco) , the difference in the dis- 

Cu edge 

persion corrections at the Cu edge between the SCO and 



LCO layers. Even for the SL higher doping levels, holes 
do not appear to empty unusual Cu states, unlike the case 
of YBCO.— Small SL imperfections shift the peaks from 
integer values. The line shape changes with L from a dip 
on resonance (low L) to a peak (high L). This change is 
not due to absorption (the calculated absorption depth 
is much larger than the SL thickness), refraction (from 
two-dimensional profiles, data not shown) or multiple- 
scattering at low L (similar effects are seen at high L in 
manganite SL, data not shown). 



D. Scattering at the O edge 

RSXS measurements at the O edge (Figs. 4 and 5) 
probe the spatial distribution of holes in the SL. Con- 
sistent with FY measurements (Sec. II B), the absence 
of a scattering peak at UHB shows that the UHB states 
have been uniformly removed in all layers. The scatter- 
ing contrast in Fig. 4 occurs mostly at the two energies 
(MCP and SHP), where SL-A and SL-C differed in the 
FY (Fig. 2). Because scattering peaks at SL reflections 
(Fig. 5), the charge density is modulated at MCP and 
SHP energies between SCO and LCO layers of the SL. 

Sharpening of the features is observed at higher L (Fig. 
4). The MCP and SHP peaks are well-separated (the ad- 
ditional splitting of the MCP in SL-C will be discussed 
in Sec. HI B). It is difficult to explain this clear sepa- 
ration with a difference in the energy of the same MCP 
hole state in the LCO and SCO layers, given the inherent 
interface roughness of the structure. This is strong ev- 
idence for considering SHP a qualitatively distinct hole 
state. 

The low L line shapes for SL-A, SL-B (not shown) 
and SL-C are more similar than the high L line shapes 
because of the reduced importance of roughness (Sec. HI 
B). At low L, the energy profiles show increased intensity 
between MCP and SHP [Figs. 4 and 5(a)], unlike the 
La M5 and M4 edges (Fig. 3). Ahhough the O edge 
line shapes are more complex, they can be qualitatively 
analyzed with the same model as for the La and Cu edges 
(Sec. HI A). The difference in the MCP and SHP hue 
shapes at low L is due to different interference conditions 
with the surface reflection (Sec. HI B). 

Vacancies are presenlji^ in bulk SCO as well as at cer- 
tain interfaces, even when the bulk materials do not con- 
tain vacancies4i If there were vacancies in the SCO lay- 
ers or at interfaces, these would give a peak in scatter- 
ing at all O edge energies. However, we do not observe 
large features above the SHP energy that would indi- 
cate a structural difference (either vacancies or oxygen 
atoms at interstitial sites), between the LCO and SCO 
layers (Fig. 5). Since the samples were annealed in 
ozone, the SCO layer thickness is relatively small, and 
there is little SL scattering at the O edge other than 
at MCP and SHP (Fig. 5), we will neglect vacancies 
in the following analysis, that is (5 « 0. The number 
of oxygen atoms is approximately the same throughout 
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FIG. 4: (Color online) Line shapes at constant L for differ- 
ent SL peaks (solid line), compared to substrate reflectivity 
(dotted line) at the O edge, for SL-A (a-c) and SL-C (d-g). 
Results on SL-B (not shown) were similar to those for SL-A, 
supporting the observation that small structural irregularities 
do not affect our conclusions. The scans have been normalized 
to unity below the edge. SL-A and SL-C have different scat- 
tering line shapes at high L because of different SL roughness 
and structure (Sec. Ill B). Measurements on the substrate, 
approximating the surface reflectivity So (Sec. Ill A), show 
that the surface reflectivity resonant contribution gets neg- 
ligible at higher L, where the substrate reflectivity becomes 
more featureless. Vertical dashed lines are guides to the eye. 



the SL, while their valence is different in the SCO and 
LCO layers, giving the scattering contrast. Therefore, 
SL scattering at the MCP and SHP energies probes the 
difference between SCO and LCO layers form factors 

jj.(MCP,SHP) //■ J. \ 

^JsL = [jsco - jLCo) , measuring 

w=MCP,SHP 

the change in the dispersion corrections with doping. 



III. DISCUSSION 
A. Scattering model 

To analyze the line shapes at the O edge, we develop a 
model for the interference between the surface and a SL 
reflection. The X-ray scattering intensity is / = ^l^p , 
with the structure factor given by: 



^(c^,g) = ^/„(^,Q)tre^«^ 



(1) 



where io is the incident X-ray energy, Q = (0, 0, 2tiL/ csl) 
is the scattering momentum in the reflectivity geometry, 
and t" is the occupation factor in atomic plane I of ele- 
ment and valence n.^^ The form factor f{uj,Q) for soft 
X-ray momenta is f{uj,Q) = /"(Q) + /'(w) + if"{uj) « 
/°(0) + f'{^) + if'i^)- For only one source of scatterers 
(n ~ 1), the variables w and Q in the structure factor 
S{uj,Q) are separable, that is S{uj,Q) = g{uj)h{Q). This 
implies identical line shapes (up to an overall scaling fac- 
tor) for the same edge at all L. This is not observed 
at the Cu edge or for the MCP state [Figs. 3(a) and 
4]. Therefore, at least two scattering sources, interfer- 
ing in the total structure factor S{uj,L), need to be con- 
sidered. "Stray light" only contributes an overall back- 
ground level. To interfere, the contributions to S{ll!,L) 
must be coherent and of the same energy. 

With the above approximation for f{u:,Q), the struc- 
ture factor for more than one type of scatterer (the con- 
dition of a SL) becomes: 



^^■'■(.;,L)^^/^(c.)5](i: 
I 



n\ 2T^iLzi/csL 



(2) 



where p^^L) = (i")e^'^'^^'/'^''^ is the Fourier trans- 
form of the distribution {(t")} of in-plane averages of the 
occupation factors for layer I and element and valence n. 
Considering only two elements, n = A and n = B as an 
example, and {tf) + (if) = 1 for all Z, we obtain a to- 
tal structure factor S^^{uj,L) 
terms: 



-^0 



-|- S^gj^ made of two 



(3) 



and 



,(B) 



PSL ^ 



(A) 
PsL ^ 



-SJs'lPsl (4) 
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FIG. 5: (Color online) Two-dimensional plots in variables en- 
ergy E and momentum L for two selected SL peaks, (a) Scat- 
tering two-dimensional resonance profile (logarithmic scale) 
for SL-C near L=2 at the O edge shows the mirror sym- 
metry of the MCP and SHP scattering line shapes. Same 
results were obtained for SL-A and SL-B. (b) Scattering two- 
dimensional resonance profile (logarithmic scale) for SL-C 
near L=4 at the O edge shows that the MCP and SHP holes 
are qualitatively different. 



where Na, Nb, N = Na + N-q are the number of ML in 
one A layer, one B layer and one superperiod, /q^j^ is the 

average of the SL layers form factors, (5/^^ = /g — 
is the contrast between the form factors of the SL layers, 
and the substrate contribution has been neglected. 

A complex oxide SL has more than two constitutive 
elements. We neglect for soft X-ray momenta the small 
difference in zi between the SrO/LaO and Cu02 planes, 
or equivalently between oxygen sites 0(1) (in-plane) and 
0(2) (apical), within 1 ML. In this case, /o becomes the 
average total form factor and Eqs. 3-4 can be applied by 
replacing Pg^'^^ with 7Va,b with A^lco.sco, 

and /a,b with /lco.SCO) the total form factors of 1 ML. 

The indices i and j are dotted with the light polariza- 



tion vectors and e^niUai as / = e*j,„^if^ej,iniuai- 

Only the in-plane /^^ remains at the Cu edge, at MCP, as 
well as at SHP from the similar angular dependence of FY 
(Sec. II B). In this case, the angular dependence given by 
the double product simplifies to cx f* f i^^^iP^'ej^initiai (x 
/^^sin^(6'), with an additional factor [sin^(6')] that can be 
absorbed into the arbitrary units. 

The momentum dependence of 5*0 and Ssl in Eqs. 3- 

4 is contained in the functions pq^sl, P^sl'~'^ ^^'^ Ps^l^*^^- 
For a SCO-LCO superlattice with no roughness, they are 
(r is the number of repeats): 



Pq,SL = > e 

LAU 



2-jriLzi/csL _ 



1 - e 



2niLr 



„(LCO) 



I _ ^2-KiL/N 



(5) 



l,LCO 



and 



^(SCO) 
PSL 



l,SCO 



PO,SL 



^(LCO) 
PSL 



(7) 



Except near L — mN, where m is an integer, the func- 
tions p^g^'^^ and p^s^'^^ are related as Pg'^*^' 
Therefore, they are out-of-phase: 



^(SCO) 
PSL 



Arg[p('^^^°)(L)] 



Arg[p(f°)(L)] 



(8) 



where Arg[p^'^^^°)(L)] = ^L(Nlco - 1)/N. 

The Sq term in Eq. 3 gives the reflection from the 
discontinuity in the index of refraction n, or equivalently 
in the average form factor fo,SL, at the sample surface. 
For a SCO-LCO SL, fo,sL = %^(2/La) + ^(2/sr) + 
fcu + 4/o , which can be calculated from tabulated^ val- 
ues. /o,SL is made mostly of non-resonant terms at the 
Cu and O edges (but not at the La edge^) and, because 
of this, does not depend strongly on energy or polariza- 
tion. Considering the substrate extends the sum in Eq. 
2 to an infinite number of layers, with the function po,SL 
replaced by poiL) = 1/(1 - e^^^*^/^) and fo,SL replaced 
by /o = (1 - T)fo,sL + rfo.LSAO, where fo.LSAO is the 
substrate contribution and t(6, 26, w) is a weighting fac- 
tor which depends on the scattering geometry {0, 20) and 
energy (cu) through a variable absorption depth. 

The Ssl term is the reflection from the SL modulation. 
The modulation of the La and Sr numbers in the SL also 
gives modulations of the form factors at the Cu (Sf^^) 

and O ((J/s^^^'^^^"*) edges. Therefore, the interference 
at different edges will be analyzed in terms of two main 
components: a SL reflection Ssl and a surface reflection 
Sq [Fig. 6(a)]. Fig. 6(b) shows the functions |po(-^)| and 
\psl{L)\ for a 8 X [2 X SCO 7 x LCO] structure. The 
L dependence of \po\ and \psl\ is very different; while 
I Pol has peaks at L = mN, where m is an integer, \psl\ 
follows a momentum dependence with pronounced peaks 
at integer L ^ mN. \psl\ is therefore much higher than 
IpoI at these L values: \psl\ >> \po\- It would appear 
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FIG. 6: (Color online) (a) Sketch of the two components in 
the structure factor S: the surface So, from the discontinuity 
between the vacuum and the surface, and Ssl , from the mod- 
ulation SfsL between SCO and LCO layers. The circles show 
the location of the MCP and SHP holes doped in the LCO and 
SCO layers, if the SHP is scattering from a hole on an orbital 
at the apical 0(2) sites (Sec. II B). A few examples of atomic 
planes in the SCO and LCO layers are indicated, (b) Calcu- 
lations of \pslAL)\ and \psL,i{L)\ for a [2 x SCO 7 x LCO] 
SL with no roughness and 8 and 1 repeats, compared to the 
surface |po(i)|. 



that the surface reflection can be neglected near SL reflec- 
tions. However, |/o| is considerably larger than \SfsL{^)\ 
at certain edges, with the SL and surface structure fac- 
tors comparable in magnitude \Ssl\ ^ \So\- 

Including both terms, the intensity of X-ray scattering 



I = A\psl\ 



SL 



foPo 



PSL 



^^\5fsL{u:) + ZL\ (9) 



with IpslP (independent of lo) absorbed into the arbi- 
trary units. It is not possible at present to calculate from 
first principles the energy- dependent form factor 5/sl(w) 
for a correlated oxide SL. However, FY measurements 
(Fig. 2) show that the energies of MCP and SHP do not 
depend on the environment of the oxygen atom, which 
suggests a simpUfied model. Therefore, the line shapes at 
different edges will be modeled with a harmonic oscillator 
functional dependence. The function used for the imag- 
inary part of the form factor is Sfgi^{uj) = j^zr^^rff^- 
This is Kramers-Kronig transformed to obtain the real 
part SfsLi^) [Fig. 7(a)]. Fig. 7(b) shows |<5/sl + zl| for 
different z^. 

The interference between the surface Sq and superlat- 
tice Ssl terms can qualitatively describe the scattering 
line shapes at all edges. At the La edge, {Sfg'j^^l is com- 
parable to |/o|, I -Ssl I >> |S'o| and little change of the 
resonance line shapes with L is expected. Indeed, there 
is little change with L at the La edge [Fig. 3(a)], except 
at low L. Measurements on the substrate at the La edge 
[Fig. 3(a)] indicate a strong resonance in /o.lsao, and 
therefore /o, which is responsible for the variation in the 
line shape at low L. Refraction effects shift the peaks and 
complicate the interpretation of the measurements at low 
L; because of this, the measurements at the La edge will 
be discussed in detail separately The scattering line 
shape at the M5 and M4 edges has the same shape be- 
cause it originates in states on the same La atoms; this 
is in contrast to the scattering at the O edge at low L, 
where the line shape at the MCP and SHP resonances is 
different (Sec. Ill B). 

At the Cu edge, 6f^^^^^ is the change in the disper- 
sion corrections with doping and, because of the small 
modulation in the Cu valence, ISfg^"*] « \fo\ and 
I'5'slI << I'S'ol- Strong interference effects are expected 
and observed [Fig. 3(b)]. The scattering shows a strong 
surface contribution to reflectivity. The non-resonant 
ifo,SL and /q.lsao) form factors have large imaginary 
parts /o because of the neighboring La edge and the line 
shapes resemble those of panels (2) and (8) in Fig. 7(b). 

The O edge scattering, with \Ssl\ ^ \So\, is an inter- 
mediate case. (5fg^^^'^^^-* is relatively large and /q' is 
smaller than at the Cu edge, with the line shapes in a 
different region of Fig. 7(b): for low L the MCP and 
SHP are of the type (4) and (6), for high L of the type 
(5) and (2) . Scattering line shapes are the most variable, 
with large complex changes with L. The scattering at 
the O edge is analyzed in more detail in the next section. 
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FIG. 7: (Color online) (a) Functions if^L and 5/^^ ob- 
tained by Kramers-Kronig transform. The parameters were 
a = 2.5 X 10"^ (top) and a = -2.5 x 10"^ (bottom), 
So = 528.5 eV and F = 1 eV. (b) Line shapes j5/sL -\-Zh \ cal- 
culated for different interference conditions, arranged in the 
(2^, 2^) plane according to the zl term used, show the vari- 
ety of behaviors at the Cu edge [panels (2) and (8)], La edge 
[panel (4)] and O edge [panels (4) and (6)]. 



B. Locus of MCP and SHP scatterers 

Interference between the surface and the SL reflections 
determines the line shape. This can be used to find the 
locus of the MCP and SHP scatterers within the SL struc- 
ture. 

As pointed out, the functions cannot 
be calculated. However, the change with doping in 
the dispersion corrections at MCP and SHP energies 
can be estimated by looking at the difference in FY 
{SfpY — f'sL~A ~ f'sL-c) between two differently-doped 
SL, which gives the scattering contrast between SL layers 
(SfsL = f'sco - f'lco) in one SL (either SL-A or SL-C). 
Specifically, when comparing SL-A with SL-C, FY inten- 
sity increases concurrently at MCP and SHP [Fig. 2(a)]. 
Therefore, for both MCP and SHP energies, the differ- 
ence in the dispersion corrections in the SCO and LCO 



layers is positive: 

(5/^'i(MCP, SHP) oc (5/^y(MCP, SHP) > (10) 

The Kramers-Kronig transform of a Sfgj^ > peak gives 
the shape shown in Fig. 7(a) (top panel) for \SfsL\ = 

However, measurements show line shapes for the MCP 
and SHP states with a mirror symmetry at L — 1 and 
L = 2 (Figs. 4 and 5). It is necessary therefore to con- 
sider the contribution to scattering of z^, the other factor 
in Eq. 9. Measured line shapes for MCP [of type (3)-(6)- 
(9), right side] and SHP [of type (l)-(4)-(7), left side] are 
on opposite sides of the complex {z'^,z'[) plane in Fig. 
7(b). Therefore, ^['^cp.shp) ^^^p gjjp ^^^^ 
posite phases. 

Considering the factors combined in 2;^^^^^'^^^^^ since 
po is the same at MCP and SHP, and /o cannot change 
sign over the 4 eV between the MCP to SHP ener- 
gies, as evidenced by the weak energy dependence of the 
substrate reflection (Fig. 4) and the large non-resonant 
component of /o, the difference between 2:^'^'-'^'^^^' must 
be in pg'^^ , the Fourier transforms of the dis- 
tributions {(i/^^^'')} and {(tp^^"*)} of the occupation 
factors at the two energies. Because and psL are in- 
versely proportional, psL also has opposite phases at the 
MCP and SHP energies: 



Arg[p^^^^\L)] 



Arg[p(fP)(L)] 



(11) 



Therefore, the distributions pg'^^^'^^^''(L) are out-of- 
phase, with the spatial distributions of the MCP and 
SHP holes (Eq. 11) related in the same way as the spa- 
tial distributions of the LCO and SCO layers in the SL 
(Eq. 8). Then, it is necessary that the MCP and SHP 
hole distributions peak in different layers. 

This conclusion, using measurements at low L, where 
roughness effects are less important, is independent of 
the thickness of LCO and SCO layers in one superperiod 
or roughness amplitude. Roughness effects become more 
important at higher L, where the MCP and SHP line 
shapes resemble line shapes closer to the origin of the 
complex plane in Fig. 7(b). We consider the L depen- 
dence of the intensity in Eq. 9, contained in the phase 
and amplitude. The zl phases for the MCP and SHP 
holes follow those of the LCO and SCO layers, which are 
linear functions of L: Arg[zL'"'~^*~'''] = Arg[fo] -fTrLNgco/N 
and Arg[z[^^^°^] = Arg[z^^°] + tt. The other contribut- 
ing factor to the change in the line shape with L, is 
the variation in the z^ amplitude. For Gaussian rough- 
ness, the ratio \po/psl\ in \zl\ depends on the sur- 
face Gs and the interface di roughness as [po/psil = 
|po/psL|ideaii?(Qz), where R{Q,) = e'^-'— and 
Qz = iTvL/csL- The factor R{Qz) becomes increasingly 
important at higher L. Therefore, the linear increase of 
the zl phase with L and the reduction in the amplitude 
of [2;i[ = \zL\\dcaiR{Qz) at higher L for ag > cti combine 
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to make the vector zl spiral inward in the (z^, z'[) plane 
of Fig. 7(b) with increasing L, systematically sampling 
different line shapes at the same edge. 

Fig. 4(c) for SL-A, with the first peak in the split 
MCP almost absent, shows the interference line shape in 
a slightly different location in the plane of Fig. 7(b) than 
the location corresponding to the line shape in Fig. 4(g) 
for SL-C, where both peaks are visible at MCP. This di- 
vergence between measurements for SL-A and SL-C with 
increasing L (Fig. 4) can now be explained by consider- 
ing the different roughness of SL-A and SL-C. SL-A has a 
larger surface ct^ than SL-C (Sec. II A) but similar inter- 
face roughness Ci. Therefore, R{Q) is smaller for SL-A, 
which makes the end of the zl vectors for SL-A and SL-C 
follow slightly different trajectories with increasing L in 
the plane of Fig. 7(b). 

The interference between Ssl and 5*0, resembling 
multi-wavelength anomalous diffraction in molecular 
crystallography, allows determining the superlattice 
MCP and SHP hole distributions on a relatively large 
scale. To confirm the FY measurements of Sec. II B and 
the site of the SHP state within 1 ML, a fit with the inter- 
ference model of Sec. Ill A would be necessary for higher 
wave vectors (L > 2), giving better spatial resolution. In 
future work, fitting the evolution of measured interfer- 
ence line shapes (Figs. 4 and 5) with L would also allow 
obtaining (5/s^^^'^^^\ Us and the energy-dependent Ci 
(on which the zl amplitude depends), and A^sco or Nj^co 
(on which the z^ phase depends). 

The MCP and SHP hole distributions, maximized in 
different layers, are consistent with the observation that 
the MCP holes are mobile^^ and that remotely doping a 
new distinct type of holes (SHP) in the LCO layers from 
the Sr in the SCO layers is an unlikely strong long-range 
process. In contrast, the SHP holes remain centered on 
the SCO layers. The relatively smaller number of MCP 
holes in the SCO layers is consistent with the saturation 
of FY in Ref. 25 and observations in Ref. 32, where 
the apical hole site is favored at high doping. The dis- 
tributions of MCP and SHP holes are illustrated in Fig. 
6(a). 

IV. CONCLUSION 

Using X-ray absorption and resonant soft X-ray scat- 
tering, we found that doping in SCO-LCO superlattices 



empties two distinct oxygen hole states. The location of 
sub-surface scatterers within the superlattice was deter- 
mined using the interference with the surface reflection, 
with the distributions of the two hole states maximized 
in different layers. 

X-ray absorption measurements suggest that the hole 
state at the higher energy is on apical oxygen atoms and 
polarized in the a — b plane. Since the creation of a 
vacancy removes two holes, this would suggest that the 
vacancies in bulk SCO are at the apical sites as well. 

The increase in the density of the unoccupied states at 
SHP has to come from an occupied state. Photoemission 
experiments to determine if a change with doping occurs 
in the density of states of occupied states that mirrors 
the change with doping seen at SHP in this study would 
be interesting. 

Calculations showed that making in-plane axial or- 
bitals more localized, by gradually spatially removing the 
apical oxygen atoms, increases Tc,^ and that, conversely, 
unoccupied apical oxygen orbitals lower Tc. If similar 
ideas are applied to the SHP state, the doping of SHP 
holes might be responsible for the decrease of Tc and the 
hole pair breaking at a; > XopHmat^, and their partial 
removal, with the appearance of vacancies in SCO near 
X = 2, for the high of bulk SCO. 
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